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Laser filamentation is responsible for the deposition of a significant part of the laser pulse energy 
in the propagation medium. We found that using terawatt laser pulses and relatively tight focusing 
conditions in air, resulting in a bundle of co-propagating multifilaments, more than 60 % of the 
pulses energy is transferred to the medium, eventually degrading into heat. This results in a strong 
hydrodynamic reaction of air with the generation of shock waves and associated underdense channels 
for each short-scale filament. In the focal zone, where filaments are close to each other, these discrete 
channels eventually merge to form a single cylindrical low-density tube over a ~ 1 p.s timescale. We 
measured the maximum lineic deposited energy to be more than 1 J • m“^. 


PACS numbers: 42.65.Jx, 47.35.Rs, 51.20.-|-d 

I. INTRODUCTION 

Laser filamentation is a propagation regime for ultra- 
short laser pulses in transparent media reached when 
their peak power exceeds a critical power Per, which is 
about 5 GW in air at a wavelength of 800 nm. It results 
from a complex interplay between the optical Kerr effect, 
nonlinear energy absorption due to multiphoton and tun¬ 
nel ionization and plasma defocusing [T]. While propa¬ 
gating in this regime, the laser pulse is able to maintain 
a high intensity over several Rayleigh lengths. Filamen¬ 
tation can deposit a significant part of the pulse energy 
in the propagation medium. Plasma generation itself, to¬ 
gether with plasma absorption, form the main channel for 
energy deposition. Stimulated rotational Raman scatter¬ 
ing on air molecules can also account for energy depletion 
of the laser pulse mil- This energy, initially stored as po¬ 
tential and kinetic energy of plasma free electrons and as 
rotational energy of air molecules, is eventually converted 
into gas thermal energy after plasma recombination and 
rotational thermalization, occurring over ~ 1 ns [5] and 
^ 100 ps timescales [5], respectively. The resulting heat¬ 
ing can range from ~ 100 K 0 to more than 1000 K 
1^, depending on experimental conditions. As air ther¬ 
mal conductivity is low, the system relaxes by launching 
an outward-propagating cylindrical pressure and density 
wave, expelling matter from the center of the channel 
and leaving a low-density tube. This underdense chan¬ 
nel then slowly evolves by thermal diffusion, which can 
take times up to several tens of milliseconds [Rj. Such 
laser-induced air hydrodynamics are particularly inter¬ 
esting because they pave the way for the development of 
remotely generated, long-lived virtual optical structures 
such as optical waveguides 0110]. They could also be 
used in the prospect of enhanced aircraft aerodynamics 
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m- As the efficiency of these structures is strongly de¬ 
pendent on the amplitude of laser-induced hydrodynam¬ 
ics effects, that is on local deposited energy, it is crucial 
to be able to increase the density of deposited energy at 
the required level. 

In this Article, we first study the influence of laser 
parameters on energy deposition. We found that us¬ 
ing the shortest and most energetic pulse with a mod¬ 
erately strong focusing yields the best results overall, 
with more than 60 % of the laser pulse energy being 
converted into heat at the terawatt peak power level. In 
this case the laser pulse peak power is well above Per, 
and the beam breaks down to form many co-propagating 
filaments due to modulational instabilities [T^ [T^ . This 
effect results in the pulse energy being deposited in com¬ 
peting short-scale filaments instead of being concentrated 
in a single structure, leading to an overall degradation 
of local energy deposition. Investigating the air hydro¬ 
dynamics response by means of transverse interferom¬ 
etry, we witnessed a spontaneous evolution of the sys¬ 
tem from discrete underdense channels to a cylindrically- 
symmetric form with a single, large channel and an asso¬ 
ciated unique shock wave, partially negating the delete¬ 
rious effect of multifilamentation. In optimal conditions, 
this channel lasts for more than 100 ms. Using a sono¬ 
graphic technique [HES], we estimated the peak lineic 
deposited energy to be 1.3 J ■ m“^. 


II. OPTIMIZING ENERGY DEPOSITION 

We first determined optimal experimental conditions 
to get the maximum absolute energy deposition from fil¬ 
amentation. To this purpose, we used a direct measure¬ 
ment of the laser pulse energy right after the focusing lens 
and ^10 cm after the end of the laser-induced plasma 
by means of a Joule meter (model QE50LP-H-MP from 
Gentec-EO). The difference between input and output 
energy was considered as deposited energy, which is a 
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good approximation since energy losses due to Rayleigh 
and Thomson scattering are negligible. Such measure¬ 
ments were repeatedly done while varying laser pulse en¬ 
ergy, laser pulse duration and focusing conditions. 

A. Laser pulse energy 

The first investigated laser parameter was laser pulse 
energy. For this study a moderately strong focusing at 
//30 was used with a 50 fs laser pulse at 800 nm. Evo¬ 
lution of deposited energy with input energy is displayed 
in figure [2 



input energy (mJ) 

Figure 1. (Color online) Evolution of the total (red circles) 
and relative (blue squares) deposited energy with input energy 
from the filamentation of a 50 fs laser pulse at 800 nm focused 
at //30. Error bars correspond to a confidence interval of one 
standard deviation evaluated over 500 shots. 

In this figure are plotted two curves: first, the total 
deposited energy (red circles), that is the difference be¬ 
tween input and output energies and second, the relative 
deposited energy (blue squares), that is deposited energy 
normalized by input energy. The global trend is an in¬ 
crease of deposited energy with input energy. Deposited 
energy evolves quasi linearly at higher laser energies. As 
for relative deposited energy, it can roughly be divided 
in three parts: the first one, up to Ei^ ~ 20 mJ, is 
characterized by a very steep increase, reaching up to 
35 %. Past this point, relative deposited energy still 
increases, but much more slowly. Eventually, beyond 
Ein ~ 110 mJ, it saturates around 60 %. We attribute 
the transition observed around 20 mJ, that is a peak 
power of about 80Per, to the threshold between the full 
beam self-focusing leading to an initial single filament 
generation and the local self-focusing inside the beam re¬ 
sulting in multifilament formation before the nonlinear 
focus, as described by Fibich and co-authors [13] ■ This 
means that as soon as the beam breaks up into short- 
scale filaments before the focal zone, where the majority 
of energy deposition occurs, increasing input energy leads 
to a redistribution of this energy over many competing 
filaments, only slightly increasing deposited energy for 
each of them. Conversely, as long as only one filament 
is concerned, any increase in energy will be channeled 


through this single structure and lead to a dramatically 
stronger energy deposition. 


B. Laser pulse duration 

Influence of the laser pulse duration on energy deposi¬ 
tion was investigated by fixing input energy and focusing 
conditions and by detuning the laser compressor, impart¬ 
ing temporal chirp to the pulse. 



Figure 2. (Color online) Evolution of the relative deposited 
energy with pulse duration in the case of positive (red cir¬ 
cles) and negative chirp (blue squares). Filamentation was 
generated by a 164 mJ laser pulse at 800 nm focused at //30. 
Error bars correspond to a confidence interval of one standard 
deviation evaluated over 500 shots. 

Results displayed in figure are unambiguous: pulse 
duration is a very sensitive parameter. As soon as it devi¬ 
ates from the minimum value of 50 fs, that is from a max¬ 
imum peak power, deposited energy quickly falls, being 
reduced from more than 60 % to less than 20 % at 1 ps, 
and to only a few percents past 10 ps. It is worth men¬ 
tioning that a negative chirp yields better results than 
a positive chirp. This can be explained by the plasma 
defocusing effect affecting the back of the pulse, which is 
stronger at longer wavelengths and will therefore tend to 
decrease the intensity of positively chirped pulses more 
than for negatively chirped pulses. With a long pulse 
duration (typically > 1 ps), plasma absorption due to 
inverse Bremsstrahlung, which is also more efficient at 
longer wavelengths, can also play a role m- 

One could think that as some energy transfer channels, 
like avalanche ionization, are favored by longer pulses, it 
could positively affect energy deposition. Such behavior 
has indeed been observed during filamentation in fused 
silica ini and water (TB]. Conversely, it seems that max¬ 
imum intensity gives the best results in air, because it 
promotes a more efficient nonlinear ionization and Ra¬ 
man absorption, in the limit of pulse durations that can 
be achieved with our laser system. 
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C. Focusing conditions 

Influence of focusing conditions was also studied by 
comparing three cases: a moderate focusing at //60, a 
stronger focusing at //30 and an extreme focusing at //3, 
corresponding to the quasi-breakdown regime described 
by Kiran et al. m- Results are given in figure It 
shows that for a given input energy, energy deposition is 
maximum in the //30 case. The three cases also exhibit 
a regime transition around 20 mJ, which is more obvi¬ 
ous for weaker focusing, with a first phase where energy 
deposition quickly increases followed by a second phase 
where it rises slowly or stagnates. 



Figure 3. (Color online) Evolution of the relative deposited 
energy with input energy using a //3 (blue squares), a //30 
(red circles) or a //60 focusing (green triangles). Filamenta- 
tion was generated by a 50 fs laser pulse at 800 nm. Error bars 
correspond to a confidence interval of one standard deviation 
evaluated over 500 shots. 


transverse interferometry on laser-generated air channels 
using the same instrument presented in references m. 
[5] and 
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Interesting features observed with the strong //3 fo¬ 
cusing are that deposited energy is significantly lower 
than in the //30 case and that relative energy deposition 
reaches a maximum around 85 mJ at ^ 50 % and then 
starts to slowly decrease, going below 40 % at 165 mJ. 
This behavior might be explained by the fact that such a 
strong focusing promotes the generation of a very dense 
and short plasma with an accordingly strong defocusing 
effect, leading to a decrease of intensity and therefore to 
a less efficient energy deposition. Even though this de¬ 
crease is modest, it results in energy deposition with //3 
focusing to be lower than energy deposition in the //60 
case at high input energy. 


III. ENERGY DEPOSITION IN THE HIGH 
ENERGY LIMIT 

A. Filamentation-generated hydrodynamics 

We then concentrated on the case that yielded the most 
important energy deposition: a 165 mJ, 50 fs (3.3 TW) 
laser pulse at 800 nm focused at //30. To character¬ 
ize filamentation-induced hydrodynamics, we performed 


Figure 4. (Color online) (a): line-integrated air density pro¬ 
files at a = 1 m after the focusing lens (linear focus) for two 
different delays, displaying self-symmetrization. (b): radial 
air density profiles extracted after channel symmetrization at 
z = 1 m. (c): evolution of the shock wave radius rsw (blue 
circles) and shock wave speed usw (red squares). Filamenta- 
tion was generated by a 165 mJ, 50 fs laser pulse at 800 nm 
focused at f /SO. 

As a large number of multifilaments is generated and, 
therefore, the laser beam exhibits a strong inhomogene¬ 
ity, we did not expect the system to have a cylindri¬ 
cal symmetry, preventing us to use the Abel inversion 
to retrieve air density radial profiles. Nevertheless it 
was found that in the focal zone after a given delay, 
air density spontaneously evolves from a highly disor¬ 
ganized state, reminiscent of multifilamentation, to a 
cylindrically-symmetric state, as displayed in figure 
(a). In this last figure are plotted two line-integrated 
air density profiles taken at the linear focus. At delay 
200 ns, several structures can be seen at the center of the 
profile, undoubtedly resulting from discrete short-scale 
filaments. The picture is very different at delay 500 ns 
because at this time, line-integrated density turned into 
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a much more organized state with a good symmetry. 

Since the system becomes cylindrically-symmetric af¬ 
ter some time, Abel inversion can be used to compute 
air density profiles. Results are displayed in figure 
(b). These profiles are characterized by an outward- 
propagating cylindrical shock wave leaving a central un- 
derdense channel, much in the same way as single fil- 
amentation. In the present case however, the shock is 
broader and its amplitude is higher than in the case of 
high-energy monofilaments (see results given in reference 
[S] for comparison), while the underdense channel has ap¬ 
proximately the same depth at about 30 % of normal air 
density, but is significantly larger, reaching a full width at 
half maximum of 2 mm after 10 ps. We still recorded the 
presence of the density hole after 100 ms, the maximum 
time we could reach due to experimental constraints. By 
this time, it enlarged so much that it affected the whole 
field of view of the camera, preventing us from extract¬ 
ing air density. Shock speed usw was evaluated using 
shock radial position rsw and is plotted in figure |^(c). 
It exhibits a steep initial increase up to ~ 500 m ■ s“^ 
during the first microsecond, followed by a slower, gen¬ 
tle decrease at subsequent times, still propagating with 
a supersonic speed after 10 ps. 



y (mm) 

Figure 5. Line-integrated air density profiles recorded at 
2 = 0.96 m after the focusing lens and for different times. 
Filamentation was generated by a 165 mj, 50 fs laser pulse at 
800 nm focused at //30. 

The symmetrization process is highly dependent on the 
position along the filament bundle and, therefore, on the 
local profile of energy deposition. For instance, looking 
at the time evolution of line-integrated air density at po¬ 
sition z = 96 cm (figure [^, we see that in this case, sym¬ 
metrization takes a much longer time than at z = 1 m 
and remains incomplete even after 10 ps, but still occurs 
by 100 ps. This gives a clue about the origin of sym¬ 
metrization. We explain it as follows: multifilamentation 
results in an inhomogeneous energy deposition, which oc¬ 
curs in several separated channels. Each filament forms 
a hydrodynamic wave and the associated air underdense 
channel. If filaments are initially close enough one from 
the other, as it is the case in the focal zone, individual 
shock waves quickly interfere constructively, leading to 


the formation of a single, large shock wave. The time 
interval during which usw increases (figure|^(c)) can be 
seen as a sign of this interference process during which 
the fastest components from the strongest filaments catch 
up with the slowest components and eventually take the 
lead. Since in the focal zone at this energy level and 
with these focusing conditions the laser pulse propagates 
in the superfilamentation regime |21j . plasma and there¬ 
fore energy deposition profile are distributed in a much 
more symmetric fashion than in standard short-scale fil¬ 
amentation areas. This explains why the shock exhibits 
an almost perfect cylindrical symmetry. As this shock is 
responsible for the formation of the central underdense 
channel by ejecting matter from the center, this chan¬ 
nel shares the symmetry of the shock wave. When sym¬ 
metrization is only partial once the shock has left (as dis¬ 
played in figure]^, thermal diffusion can act over longer 
timescales to smooth out any remaining sharp density 
features, yielding a diffusive symmetrization, unlike the 
pure acoustic symmetrization process occurring near the 
linear focus. 


B. Lineic deposited energy 

Still investigating the case of a 165 mJ, 50 fs laser pulse 
focused at //30, we estimated the lineic deposited energy 
using a sonographic technique. It consists in scanning the 
acoustic emission from the filament bundle along the laser 
propagation direction by means of a microphone (model 
4138 from Briiel & Kjaer) (HKTS]. Indeed, since heating 
due to filamentation occurs much more quickly than the 
characteristic development time of air hydrodynamic re¬ 
sponse, we can consider it as an isochoric process. If we 
write the deposited energy as AC/ then the first law of 
thermodynamics reads: 

AU = CyUo [ {T{r’) - Tair) d^r^ 

= ^[ iP{^) - Pair) d^r', (1) 

where Cy is the isochoric heat capacity of an air molecule 
and no = 2.47 x 10^^ m“^ is the air density at pair = 
1.013 X 10® Pa and Tair = 300 K. In the case of an initial 
pressure profile with cylindrical symmetry, this equation 
can be simplified as: 

AU « ^ [ {Pmax{z) -Pa*r)?’o(-z)^ dz, (2) 
Kb Jr 

where ro{z) is the half width at half maximum (HWHM) 
of the temperature profile at z. Assuming that vg does 
not vary significantly along the acoustic source, we find: 

AU (X {p 

max {z) - Pair) dz. (3) 

JR 

Recording the z-evolution of the peak pressure with 
a microphone in the case of an initially cylindrically- 
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symmetric pressure profile then yields the lineic de¬ 
posited energy. If this method works well in the case 
of a single filament [22], it cannot be directly applied in 
the present case because this requirement is not fulfilled. 
The symmetrization process in the focal zone can be used 
at our advantage. As the system spontaneously evolve to 
this cylindrically-symmetric state, we can consider it as if 
it started from a cylindrical pressure profile with a similar 
initial energy. We then make the assumption that most 
of energy deposition occurs around the focal zone where 
this effect is at its strongest, enabling us to retrieve the 
lineic deposited energy following equation 


tegral of the corresponding curve with deposited energy 
following equation ^ . Results are plotted in figure 
(b). They exhibit a zone about 15 cm long where lineic 
deposited energy is above 1 mJ ■ cm“^, reaching a peak 
value of 13 mJ ■ cm“^ shortly before the linear focus, 
confirming that most of the energy deposition takes place 
in the focal zone. As a comparison, high-energy single 
filamentation presented in reference jS] only reached a 
400 pJ ■ cm“^ maximum energy deposition. 

C. Physical investigation of energy deposition 


(a) 



t(|JS) 


(b) 



Figure 6. (a): example of microphone signal recorded at 

a = 80 cm, displaying only a single acoustic wave, (b): spatial 
evolution of the lineic deposited energy from a 165 mJ, 50 fs 
laser pulse at 800 nm focused at //30. This curve was ob¬ 
tained from sonographic measurements of the bundle acoustic 
emission and by equating the z-integral of the corresponding 
scan with total deposited energy. Error bars correspond to a 
confidence interval of one standard deviation evaluated over 
200 shots. 

As seen in figure]^ (a), the filament bundle already gen¬ 
erates a single acoustic wave as far as 20 cm away from 
the linear focus, showing that at least a partial acoustic 
symmetrization occurs at this position and strengthen¬ 
ing our hypothesis. We therefore performed a full acous¬ 
tic scan along the filament bundle and equated the in- 


investigating the underlying physics of laser energy de¬ 
position asks for a good knowledge of the different vectors 
for energy absorption. Plasma generation is the easiest 
one to characterize. To this purpose we used a spectro¬ 
scopic analysis of the plasma luminescence. It was shown 
previously by Theberge et al. that emission from the first 
negative system of the NJ cation can be used as a plasma 
diagnostic following: 

[ L{z) dz cx W, (4) 

Jm. 

where L is the luminescence signal at position z along 
the multifilament bundle and Ng the total number of free 
electrons in the plasma [22] . 



Figure 7. (Color online) Plasma luminescence from the first 
negative system of at 391 nm (blue squares) and lineic 
energy deposition profile (red circles) along the multifilament 
bundle. Error bars correspond to a confidence interval of one 
standard deviation evaluated over 200 shots. 

This means that plasma luminescence from the first 
negative system of is proportional to the area- 
integrated electron density: 

L(z) oc / ne{r,6,z)r drdO. (5) 

Jm? 

We recorded plasma luminescence along the multifila¬ 
ment bundle at wavelength 391 nm using a monochro¬ 
mator (model H-20 from Jobin Yvon) with a ^ 1.5 nm 
spectral resolution. Results are plotted in figure [^ com¬ 
pared with the lineic deposited energy profile. 



























































6 


From this figure, one can see that the total plasma 
length is in excess of 15 cm, starting slightly before z = 
90 cm and ending between z = 105 and z = 110 cm. It is 
interesting to note that energy deposition clearly starts 
increasing before plasma luminescence, indicating that 
in the area before the focus ionization alone is not the 
main channel for energy deposition. We attribute this be¬ 
havior to the occurrence of stimulated rotational Raman 
scattering on air molecules [5]. Conversely, once in the 
focal zone, lineic deposited energy and plasma lumines¬ 
cence behave quite similarly, decreasing simultaneously 
and almost completely disappearing past z = 105 cm. 
This good agreement shows ionization becomes the main 
channel for energy deposition close to the focus. 

Assuming that the peak energy deposition is reached 
mostly due to ionization of O 2 molecules [T], the peak 
lineic deposited energy is given by: 


dAU 

dz 


= / ne,max{.r,Zmax)Uo2r ^rdQ 

~ ^e,maa:('2max)Co2'^^0’ (^) 


where C /02 ~ 12 eV is the ionization potential of O 2 , 
Zmax the position along the multifilament bundle where 
the peak lineic energy deposition is reached and tq the 
HWHM of the electron density profile. Peak electron 
density in the case of single filaments generated using 
a ~ //30 focusing has been recorded [53] and numeri¬ 
cally estimated [53] to be on the order of 3 x 10^^ m“^. 
Since superfilamentation typically results in a peak elec¬ 
tron density one order of magnitude higher than in equiv¬ 
alent single filaments, with an equivalent radius in the 
focal zone of ~ 300 pm [5T], the peak lineic deposited 
energy is therefore estimated at: 


~ 15 mJ ■ cm (7) 


dAU 

dz 


which is indeed close to the measured peak value. 

Energy deposition then seems mainly limited by two 
different factors: first, electron density that can be 
reached during pulse propagation and, second, the ef¬ 
fective length of the plasma column. Given the results 
presented in this Article, it is clear that increasing the 
focusing leads to a rise of peak electron density, com¬ 
pounded by a slower shrinking of filamentation length, 
leading to an overall increase of energy deposition. How¬ 
ever a too strong focusing is detrimental to energy depo¬ 
sition because of the generation of a very localized and 
strong plasma, preventing the pulse energy from being 
efficiently focused and deposited in air. 


IV. CONCLUSION 

In this Article, we presented results from the study 
of energy deposition in air in the high energy, multifila- 
mentation regime. Influence of experimental parameters 
on total energy deposition was first investigated. It was 
found that deposited energy increases with the laser pulse 
peak power, both in terms of energy and duration. Influ¬ 
ence of pulse duration is important, with slight deviations 
from the minimum duration yielding a significantly lower 
energy deposition. Influence of pulse energy depends on 
whether the beam forms a single or multifilaments. In the 
first case, energy deposition quickly increases with input 
energy as more and more energy is channeled through the 
single structure. In the second case, deposited energy still 
increases with input energy, but does so at a far lower 
rate as excess energy is redistributed over many com¬ 
peting filaments. As for focusing conditions, we found 
that strong focusing yielded a significantly higher energy 
deposition than weak focusing. However increasing the 
focusing too much actually leads to a significantly lower 
energy deposition, due to the strong generated plasma, 
preventing from efficiently channeling laser energy at the 
focal point. Optimal focusing was evaluated at //30 for a 
165 mJ pulse energy and a 50 fs pulse duration, yielding 
the deposition of more than 60 % of the input energy. 

Filamentation-induced hydrodynamics were then stud¬ 
ied by means of transverse interferometry. As the beam 
propagates in the multifilamentation regime, initial en¬ 
ergy deposition takes place in many parallel structures 
in an inhomogeneous fashion. However in the focal zone, 
where superfilaments occur |21j . air structures sponta¬ 
neously evolve towards a quasi cylindrically symmetrical 
state after a delay typically on the order of 1 ps. We ex¬ 
plain this phenomenon by the constructive interference 
between individual shock waves generated by each su¬ 
perfilament, and the relatively cylindrical distribution of 
these superfilaments. Going away from the focal zone, 
symmetrization due to the shock wave process is only par¬ 
tial but still happens over longer timescales 100 ps) 
due to thermal diffusion. This symmetrization process 
enabled us to use a sonographic analysis to estimate the 
spatial evolution of the lineic deposited energy. It was 
found that energy densities above 1 mJ ■ cm“^ could be 
reach over 15 cm, with a peak lineic energy of more than 
1 J ■ m“^. 

These results are particularly interesting in the 
prospect of filamentation-induced air virtual optical 
structures, such as virtual waveguides. Multifilamenta¬ 
tion indeed gives birth to significantly deeper and longer- 
lived underdense channel that could be used in a similar 
way as described in reference [3]. A major drawback 
of multifilamentation comes from the inhomogeneity of 
the whole beam, especially the existence of a zone where 
no symmetrization occurs. Still, some degree of control 
can be achieved over filamentation occurrence m, which 
could help shape energy deposition profile. Also, using 
a smoother and cleaner beam spatial profile would de¬ 
crease the modulational instability of the beam and allow 
to have a higher threshold for the appearance of multifil¬ 
aments m- 
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